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ABSTRACT 

Upcoming observatories will be able to detect the kinetic Sunyaev-Zel'dovitch (kSZ) effect with un- 
precendented signal-to-noise, and cross-correlations with foreground signals such as galaxy counts are 
a promising way to extract additional cosmological information. We consider how well a tomographic 
galaxy-count cross-correlation experiment, using data from WMAP, ACT and SALT, can significantly 
constrain the properties of dark energy. We include the need to model a wide range of effects, includ- 
ing those associated with complicated baryonic physics, in our analysis. We demonstrate how much 
of the cosmological information contained in the kSZ comes from larger scales than that in the galaxy 
power spectrum, and thus how use of the kSZ can help avoid difficult systematics associated with 
non- linear and scale-dependent bias at k > Ih Mpc^ 1 . 
Subject headings: 



1. INTRODUCTION 

Future microwave experiments such as Planck, the At- 
acama Cosmology Telescope (ACT) and the South Pole 
Telescope (SPT) will map the cosmic microwave back- 
ground (CMB) with unprecedented resolution and preci- 
sion. In anticipation of this new era, much attention 
has been drawn to what can be learned from studies 
of not only the so-called "primary" CMB anisotropies 
- those arising at the surface of last scattering - but also 
"secondary" anisotropies, created by structures at much 
lower redshift. 

Secondary anisotropies arise at all angular scales; the 
largest secondary anisotropy at the arcminute scales 
probed by the experiments above is the thermal Sunyaey - 
Zel'dovitch (tSZ) effect l|Sunvaev fc Zeldovichl IT970TL 
caused by scattering of CMB photons from the hot, ion- 
ized gas associated with clusters. The tSZ has a unique, 
non-thermal signature which, in principle, allows it to be 
isolated from other contributions. 

After removal of the tSZ, the dominant source of ar- 
cminute anisotropy is expected to arise from the kinetic 
Sunyaev-Zel'dovich (kSZ) effect, caused by bulk flows in 
the same ionized ga s, and the subject of this p aper. It 
was first studied bv lOstriker fc Vishniad l)1986|) and for 
this reason is also known as the Ostriker-Vishniac ef- 
fect (OV) 3 . Other important contributors to arcminute- 
scale anisotropies include the gravitational lensing of 
prim ary anisotropies an d, potentially, patchy reioniza- 
tion ijSantos et al.l20 03L e.g.); beyond £ of approximately 
3000, the amplitude of the primary anisotropies t hem- 
selves is subdominant because of Silk damping ijSilkl 
1968). 

In linear perturbation theory, the magnitude of veloc- 
ity perturbations, and thus the kSZ, is directly related to 
both the growth rate of structure and the Hubble con- 
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3 Because lOstriker &: Vishniad <198(3) assumed both the velocity 
and density fields were linear, studies often refer to the "OV" effect 
as that arising from purely linear fluctuations, reserving the term 
"kSZ" for calculations that include a non-linear treatment of the 
density field. 
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where D is the growth factor at late times. A measure- 
ment of the amplitude of the kSZ and, in particular, a 
measurement of the amplitude as a function of redshift, 
would provide a new way to probe the nature of dark 
energy, which influences both the geometry of spacetime 
and the growth-rate of structure. 

In this paper, we consider how well a measurement 
of the kSZ, in cross-correlation with a redshift-binncd 
galaxy count signal such as might be provided by the up- 
coming Southern African Large Telescope (SALT), might 
constrain properties of dark energy, including both the 
equation of state parameter w, and its derivative with 
respect to scale factor, dw/da. 

The kSZ is not a "clean" signal: it arises in part from 
the non-linear regime of structure formation, and de- 
pends on the details of baryonic physics on megaparsec 
scales. We first review, in Sees. |3] and 03 the calcula- 
tion of the angular power spectrum from the underlying 
three-dimensional density fields. 

In Sees. 0] and |3J we examine some of the approxima- 
tions made in this calculation, comparing our results with 
those from simulations that include gas dynamics, and we 
introduce additional parameters into our model to take 
into account the need to model the effects of baryonic 
physic s. This is an extension of the work of iDore et alJ 
( 2004) , who presented the first semi-analytic calculation 
of a cross-correlation between the kSZ and a density 
tracer (in that case, the weak lensing of galaxy images.) 

In Sec. El we describe the Fisher matrix analysis we use 
to determine how the combination of ACT and SALT 
data, supplemented by large angular-scale data from 
the Wilkinson Microwave Anisotropy Probe (WMAP), 
can determine cosmological parameters including w and 
dw/da. We present our results in Sec. In Sec. |H1 
we examine how the use of the kSZ effect may allow us 
to partially circumvent difficult questions of small-scale 
physics. We discuss our results, and suggest future av- 
enues for investigation, in Sec. ED 

2. THE KSZ AND GALAXY COUNT SAMPLES 
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In order to compute both the ACT kSZ signal and its 
cross-correlation with a SALT ga laxy-count s ample , we 
follow and extend the methods of lDore 

The kSZ signal is caused by the scattering of CMB 
photons off of moving ionized gas; unlike the tSZ it does 
not lead to a spectral distortion but rather only, to first 
order in v/c, a temperature increment or decrement de- 
pending on whether the halo is moving towards or away 
from the observer. The fractional change in temperature, 

e = at/To, is 
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where 7 is the unit vector pointing away from the ob- 
server and the momentum density, q, is defined as 

q(iv, »?) = [! + S(frTf)]€(rn, 77), 

<7, the visibility function, is defined as 

g{ri) = x e T H {\ + z) 2 e~ r , 

and thi the Thompson-scattering optical depth to the 
Hubble distance today, is 

t h =0.0691(1-1^6/1, 

where Yp is the primordial Helium fraction. As discussed 
in Sec.Iolbelow, we describe the reionization linear 
function of redshift, zero before some early z T i. 

For our fiducial SALT galaxy sample, we assume an 
idealized tracer with a linear bias on all scales. As dis- 
cussed below, it is certainly possible to conduct an anal- 
ysis with a less ideal tracer - and some methods of kSZ- 
galaxy cross-correlation require detailed attention to bi- 
asing - but for simplicity, and in ignorance of the actual 
biases that will be uncovered in a SALT sample, we do 
not consider this in detail. 

We take our sample to have 6 x 10 5 galaxies over a 
total area of 100 square degrees, with a median redshift of 
0.7, distributed log-normally with a first-moment equal 
to 0.3 and with sufficient redshift information to allow 
the sample to be divided for the purposes of tomography 
into ten redshift bins with equal numbers of galaxies. 

We define Ws t i(rj) as the window, or galaxy selection 
function of our sample; the subscript i refers to one of 
ten equal number bins; W is proportional to the num- 
ber density of galaxies within the bin, and zero outside. 
For conciseness, we will omit this subscript in our later 
equations. 

For the ACT experiment, we assume a #fwhm of 1.7 
arcminutes and noise cr p i X of 2.0 /iK covering the same 
sky region as the SALT sample. 

3. COMPUTING ANGULAR POWER SPECTRA 

The kSZ effect only becomes significant at very small 
angular scales (I ~ 10 3 ). We thus use the Limber 
(Kaiser) approximation when projecting the various 3- 
d density fields responsible for both the kSZ and galaxy 
angular power spectra. 

Our goal here is not to produce a "final answer" to 
the question of modelling the kSZ, but rather to demon- 
strate some useful approximations. We point out and 
model, to a first approximation, where difficult questions 
concerning gas physics arise. We demonstrate a model 
that roughly represents the amplitude and shape of the 



kSZ and its cross-correlation, and the sensitivity of these 
measurements to various cosmological parameters. One 
important question that we do not address here is the 
possibility of patchy reionization. 

The kSZ effect is symmetric: an electron overdensity 
responsible for scattering a CMB photon is as likely to 
be moving towards us as away, and the temperature per- 
tur bation is proportio nal to the velocity. We again fol- 
low IDore et all lj2004T) in correlating the square of the 
(filtered) CMB temperature with the galaxy count sig- 
nal. Filtering is necessary, since squaring in real space 
amounts to a convolution in Fourier space; the kSZ sig- 
nal will be lost under the contribution from primary 
anisotropies and system noise unless the latter compo- 
nents are filtered: 

0/0 - /(f)O(f). 

The square of the filtered temperature field can then be 
written 
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and the correlation between a galaxy-count and kSZ 
squared signal is 
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The three-point function in the integral can be written 
as the projection of a three-dimensional power spectrum 
as 
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where we define the hybrid bispectrum as 

(qj(h)qj(k 2 )S(k 3 )) = (27r) 3 B q . q . s (ki,k 2 , fa) 
x8 D (ki +k 2 + fa), 
where qx, is shorthand for 3-7. Following lBuchalter et all 
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( 2000) , we ignore radial modes (the Limber approxima- 
tion), taking k z to be much less than ljr\ and finding 

{Q f {t)Q f {f)5{l"))^{2nf5 D (l + f)f(l")f{\l 
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This gives us all we need to compute C t f , defined as 

where T, the "triangle power spectrum" , is defined as 

7^T2/(w)/(|£ + q\v)Bq iq ^s(q, -k - q, k). 

As discussed in IDore etall l|200l . this can be under- 
stood as an integral over all triangles with sides (q,—k — 

q, k), where the triangles are restricted to lie on planes of 
constant redshift. It is a way to "collapse down" the rich 
set of information contained in the three-point function. 

When considering future analyses of the kSZ signal 
dealing with real data, it seems likely that a gain in 
signal-to-noise can be had by constructing a filter to 
preferentially magnify those triangle configurations sen- 
sitive to the parameters of interest, in a fashion similar 
to analyses that have bee n done in the area o f lensing of 
the primary anisotropies ijCoorav et al.ll2000j) . Since our 
goal in this paper is to get a rough estimate of the overall 
signal-to-noise as it pertains to extraction of cosmological 
parameters, we ignore this subtlety for now. 

4. THE HYBRID BISPECTRUM 

The momentum perturbation is [1 + 5(x)]v(x). In 
Fourier space, this is 

C rl 3 k' 

q{k) = v(k) + J —v(k')5(k - k'). (2) 

We are interested in the component parallel to the line 
of sight, q"j. We will work exclusively in the small angle 
approximation, so that a position on the sky is writ- 
ten 9. The kSZ autocorrelation in the two-point Limber 
(Kaiser) approximation has been worked out in many 
places; we will work out the kSZ-kSZ-galaxy correlation 
in the three-point Limber approximation. 

Determining an analytic form for the kSZ signal is a 
series of approximations. We will refrain from substitut- 
ing in definite forms for the various two- and three-point 
functions until the next section. As disc ussed in the origi- 
nal analysis l)Ostriker fe Vishniadl 198^1 . the "trick" is to 
note that, in the Limber approximation, only q modes 
with k perpendicular to the line of sight contribute. We 
can write down the "curl" component, q p , which has ve- 
locity perpendicular to its wavenumber. The magnitude 
of the line-of-sight momentum component squared is, fi- 
nally, (l/2)q p -q p , where the factor of one- half comes from 
the averaging over sin 2 4>; q p can point anywhere in the 
plane perpendicular to the line of sight. 

Within linear theory velocities are purely gradient and 
the first term in Eq. will not contribute to the curl 
component, so that only the mode coupling term Sv is 
important in each instance of q&. Many recent analyses 
have taken the density field to be non-linear, but kept 
the assumption of line ar and gradient velocity fields ( Hu 
ffll & Fryll2002r. 

Zha ng et alJ <|2004f> . however, have raised doubts as to 
whether the assumption of a linear velocity gives the full 
amplitude of the kSZ signal. In particular, they find that, 
at scales relevant to ACT, curl modes of the non-linear 
velocity field may be of similar magnitude to the gradient 



modes estimated with the linear density field. The ve- 
locity curl mode is generated by shell crossings and the 
resultant gas physics, and is probably intractable ana- 
lytically. While still considering only the Sv term, they 
suggest a phenomenological substitution of the fully non- 
linear density field into the equation for the velocity from 
linear perturbations (Eq. Q]) Although this substitution 
is hard to justify physically, it accords very well with 
their simulations, and we follow their practic e. 

Apart from this sub stituti on, Zhan get alJ l)2004|) use 
the same reasoning as IMa fc Frvl 1)2002) when comput- 
ing the kSZ power analytically, except that in the lat- 
ter paper certain approximations are later introduced to 
make the mathematics of the two-point function more 
tractable. As we show below, when computing the three- 
point function there are ten separate terms to keep track 
of; intr oducing th e ana logous approximation to that 
made in IMa fc Frvl l|2002j) significantly simplifies the cal- 
culation while still providing a good estimate of the true 
kSZ signal at small scales. 

A potent ially complicating factor has been raised by 
iPeell l)2005|) . who noted that a naive application of the- 
ory undcrpredicts even large-scale velocity correlations 
estimated when the kSZ effect is used to determine the 
proper motions of clusters: in that case, the fact that 
clusters form preferentially in already overdense regions 
demands that biasing effects in the velocity must be 
taken into account. 

In this paper, by contrast, we consider for simplicity a 
cross-correlation of the kSZ with an "ideal" galaxy sam- 
ple with linear bias, and do not preferentially select clus- 
ters. 

We now proceed to calculate the three-dimensional 
three-point function for the kSZ-kSZ-galaxy count sig- 
nal, 

(qj(ki)q^(k 2 )5(k 3 )) = ^(q p (h) ■ q p (k 2 )S(k 3 )) . 

There will be 5 x 4/2!, i.e., ten, different products of 
two- and three-point functionals to consider, but many 
of these will turn out to be zero. Schematically, and 
assuming that the connected five-point term is zero, the 
terms will look like 

(SvSvS) = (55)(vvS) + {6v)(Sv6} + (vv)(S55) 

+ (SS)(Sv5) + (vS)(5vS), (3) 

with some terms occuring more than once, with a per- 
m utation of the ( suppressed) indicies. The Ansatz 
of lDore et alJ |2004) effectively considered only the third 
term of the series. Below, we will judge the validity of 
this approximation. Let us, to begin with, examine this 
third term to determine its complete analytic form. We 
have 

P ii 3 k r H 3 k' 

(w)(sss) = (5%-kiki) jjJL] —{v\ky(k')) 

{S{k 1 -k)5{k 2 -k')S{k 3 )). 



Substituting in the velocity power spectrum and the 
density bispectrum we find this to be equal to 




{2n) 3 5 D {k 1 + k 2 + k 3 -k- k')B S ss{ki -k,k 2 - k' , k 3 ), 
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term 



{■u I (fc)^(fc'))(5(fci - k)5(k 2 - k')S(k 3 ) 
(i/(fe)<5(fel - k))(vi(k')&(k 2 - k')5(k 3 ) 
{^(jfc^jfea - k')}(6(ki - k)vi {k')8[k 3 ) 
(8(k 2 - k')vi(k'))(6(ki - fcy (fc)<5(fc 3 ) 
(<5(fci - A-)^'(fc'))(t; i (fc)5(fc2 - k')S{kz) 
{v\k)8(k 3 )){5(ki - k)5(k 2 - k')vi (%') 
(5{^3)vi{k')){5{%i - k)8(k 2 - £'y (£) 
(S(ki - k)8{k 2 - k')}^ 1 ^^ (k')S(k 3 ) 

(S(k 2 - fc')5(fc3))(« i (fc)« i (fc') 5 (^i - k) 
(Sjki - k)8{k :i ))(v i {k)vi{k l )5(k 2 - k') 



scaling 



1 


k/k 2 


k/ki 
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Fig. 1. — The ten terms of the the kSZ-kSZ-galaxy cross- 
correlation, and the scaling of their associated geometric coeffi- 
cients after integrating over k 1 . A zero implies that the term can- 
cels after integrating over the angular direction of k. 



which, upon working out the delta functions, becomes 

d 3 k 



(27TfS D (k 1 +k 2 +k 3 ) 

P vv (k)B sss (ki - k,k 2 + fc, fc 3 ). 



[1 - (ki ■ k)(k 2 ■ k)] 



We now examine the high-fci.2,3, nonlinear behaviour. 
Because the peak of the integrand will occur at some 
small k, we can, following IMa fe Frvl (|2002f) . drop terms 
of order fc/fci.,2,3- After averaging over the k orientation, 
and projecting along the line of sight, and assuming the 
other terms in Eq. [3] are negligible, we find the hybrid 
bispectrum, B q „ q . s {ki,k 2 ,k z ), to be 



1 - ifci • fc 2 ) v 2 ms Bsss(ki, k 2 , fc 3 ). (4) 
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T his analytic result is very close to the Ansatz, made 
in IDore et al. (2004), Eq. 30, which turns out to be a 
slight underestimate of our Eq. 0] 

We have a number of other terms to consider, now, to 
see if they can be eliminated in the high-fc approximation. 
The essential "trick" is that there are geometric factors 
in front of these other terms that either cancel or are of 
order fc/fci,2,3- Fig. 1 shows the ten different cases. As 
can be seen, the term that leads to Eq. 0] is the leading- 
order term in the nonlinear regime. A more precise model 
can be constructed by retaining terms of order fc/fci,2,3 
in Eq. 01 and including the three additional terms, listed 
in Fig. 1, that do not cancel by angular symmetry. 

Having reduced the hybrid bispectrum down to a for- 
mula involv ing the density bispectruim we use the fitting 
formula of iScoccimarro fe Couchmanl 1)20011) to deter- 
mine the (non-linear) Bsss; schematically, this involves 
permuted pairs of the non-linear power spectrum, P . 

B 5 s8(ki,k 2 ,k 3 ) = 2Ff{h, k 2 )P$(h)P${k 2 ) + cyclic. 

To compute the non-linear power spectrum, we use 
the halofit code, which implements the formulae of 
Smith *etaD ((2003). We modify the code to allow for 
a w(a) dependence. Note that halofit uses an analytic 
approximation to the transfer function that does not in- 
clude the phenomenon of "baryon wiggles" ; our estimate 
of the cosmological parameters thus does not include the 
additional information contained in these features. 



The function F 2 appearing above is a complicated, 
though computable, function of the parameters esti- 
mated by halofit. As described in IDore et a D 1I200I . in- 
tegrals involving B$s$ may be significantly simplified by 
invoking the same, high-fc, approximation used to derive 
Eq. The triangle power spectrum, T, defined above 
may then be written as 

A 2 T (k,z) = ^ ma Al l (k,z)E 3 (k,z) 

where the dimensionless power spectra and A^j are 
defined as 



A 2 T 



k 3 

K pnl 



and E3 is defined as in IDo re et alJ ll2004t) based on t he 
expressions found in IScoccimarro fc Couchmanl (|2001t) . 

5. MODELING BARYONIC PHYSICS 

As hot gas condenses into stars, is influenced by super- 
nova and AGN feedback, or even as it simply responds 
to its own pressure gradients, we expect the distribution 
of the baryons responsible for the kSZ effect to diverge 
from that of t he underlying dark m atter. Current X-ray 
obser vations dV ikhlini n et al.l 12005^) and cluster simula- 
tions (Kravtsov et al. 2005) find that the baryon fraction, 
/g as , approaches the cosmological value beyond r§oo, the 
radius at which the cluster is 500 times overdense com- 
pared to cosmic, but the details of this behaviour, and 
how it depends on cluster mass, have yet to be completely 
understood. 

We are thus unable to precisely model the systemat- 
ics of baryon physics as they relate to the determina- 
tion of cosmological parameters. However, it seems likely 
that a combination of computer simulations and obser- 
vational work will gradually refine our understanding of 
baryon physics to the point where systematic errors may 
be strongly constrained. In anticipation of that era, we 
include a "toy" model of baryon physics whose free pa- 
rameters are included in the final Fisher matrix analysis. 

Our model first allows the total baryon fraction to 
evolve as a function of scale factor; i.e., we take / gas (a) 
to be Q& (1 + ca), where a is the scale factor and c is a 
free parameter to be determined by the model fit. This 
is intended to model the incorporation of baryons into 
stars. 

Since the halofit code transfer function does not in- 
clude the "baryon wiggles," a changing flj, at late red- 
shifts is indistinguishable from a changing ionization 
fraction, x e . From the combination of the WMAP 
measurement of the optical depth to reionization, and 
the detection of th e Gunn-Peterson trough at z of 6.28 
(jBecker et al.ll200l|) . we know that the reionization his- 
tory of the universe is more complex than a simple step 
function at some early redshift; the inclusion of the pa- 
rameter c goes some way to including the need to "fit 
out" for these effects when determining the cosmological 
parameters. 

While we have thus included the possibility of a time- 
dependent x e , there is also the possibility of a spatial 
dependence, a phenomenon known as p atchy reionization 
and a lready detected at low redshifts l)Oh fc Furlanettol 
2005). Patchy reionization has the potential to reduce 
the signal-to-noise of the kSZ detection; studies of its 
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Fig. 2. — A comparison of the baryon (dashed line) and dark 
matter (undashed line) power spectra at redshift zero from our 
simulations. At scales smaller than approximately 7 Mpc the two 
spectra begin to diverge; our smoothing prescription (described in 
the text) allows us to accurately predict the baryon power spectrum 
from that of the dark matter down to scales of approximately 1.5 
Mpc. 



effect on kSZ detection suggest that "worst case" models 
may reduce the signal-to-noise by a factor of two. For 
simplicity, we do not consider the effect here. 

Secondly, we allow for a "smoothing scale", cr sm , the 
width of a Gaussian convolved with the dark matter dis- 
tribution to determine the baryon distribution. In accor- 
dance with both observation and simulation, the baryons 
are "smoothed" out compared with the dark matter. The 
smoothing scale not only introduces partial degeneracies 
in our final parameters, it also reduces the amplitude of 
the kSZ signal, and thus our overall signal-to-noise. 

The inclusion of a smoothing scale fits both with sim- 
ple adiabatic gas simulations, where pressure leads to gas 
profiles less concentrated than the associated dark mat- 
ter, and with more complicated simulations that include 
cooling and star formati on, which produce evi dence com- 
patible with smoothing ijKravtsov et alJ l2005'> , including 
a gre atly reduced f^s near the cluster center. Observa- 
tions ( Vikhlini n et al.l2005[) also find similar effects com- 
patible with our model. 

To determine a realistic value of cr sm , we analyse a 
200ft. -1 Mpc hyd rodynamic s i mulat ion. The Eulerian 
TVD+PM code o flTrac fc Penl l|2004ft was run with 1024 3 
grid cells and 512 3 dark matter particles. The cell spac- 
ing of 195ft -1 kpc and particle mass of 4.96 x 10 9 ft -1 M© 
provides sufficient resolution for this work. Note that this 
code h as been used previ ously in similar studies of the SZ 
effect l|Zhang et al.l2004l e.g.). Results for a 100ft -1 Mpc 
simulation with the same number of cells and particles 
are identical in the region of overlap. 

In Fig. |21 we plot the baryon and dark matter power 
spectra from this simulation at redshift zero, when we ex- 
pect effects from baryon physics to be most pronounced. 
As can be seen, the baryon power spectrum diverges from 
the dark matter at lengths of around 7 Mpc. 

Setting the smoothing scale, er sm , equivalent to a length 
of 2.2 Mpc, we can model the baryon power spectrum to 
ten percent accuracy down to lengths around 1.5 Mpc at 
redshift zero. Below this length our naive model tends to 



underpredict structure in the baryon distribution. Our 
simulations suggest that cr sm , if taken to be a comoving 
quantity, changes little over the redshift range of interest. 

We do not expect our parameters to be the "final 
answer" to the difficult problems of modeling baryonic 
physics. However, given the current uncertainty, our two 
parameters, c and er sm , covering the possibility of redshift 
and scale dependence of the ionized baryon distribution, 
seem to be a reasonable way to study both how uncer- 
tainties about the effects of baryon physics degrade our 
ability to study the properties of dark energy, and how 
studies of the kSZ effect may be a new means, compli- 
mentary to studies in the X-ray, to probe the distribution 
of gas. 

6. DETERMINING COSMOLOGICAL PARAMETERS 

We combine the following information to determine 
our set of parameters: large angular-scale WMAP data, 
small angular-scale ACT data, and a SALT galaxy-count 
sample that overlaps the ACT strip. We consider a "to- 
mographic" analysis: i.e., in our analysis we divide the 
galaxy count sample into ten redshift bins with equal 
numbers of galaxies in each. Given the log- normal galaxy 
distribution described above in Sec. |3 we are most sen- 
sitive to fluctuations around redshift of 0.7. 

As mentioned above, we require a filter, f(l), to fil- 
ter the temperature map to avoid contamination of the 
kS Z signal from "out of band" power upon squaring. As 
in iDore et al. (2004), we use the filter that maximizes 
the kSZ signal alone, 



f(l) 



£<kSZ • 



(5) 



where the "noise" term includes both instrumental noise 
and any sources of CMB power not associated with the 
kSZ. 

To characterize the performance of an ideal experi- 
ment, we compute the Fisher matrix errors; from the 
Crameo-Rao inequality, these are the best possible error 
bars that can be achieved by the above described ex- 
periments given their noise characteristics and assuming 
Gaussianity. The Fisher matrix is 



F, 



a/3 



EE 

l X.Y 
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(6) 



where A Q refers to the set of parameters we hope to de- 
termine, and X and Y refer to the twenty-one different 
measurements, i.e., temperature-temperature, galaxy- 
galaxy (ten bins), and temperature-squared-galaxy (ten 
bins.) In the case of the temperature-temperature cor- 
relation, we use WMAP data on large scales, where the 
sky coverage beats down cosmic variance, and ACT for 
small scales. 

The covariance matrix has off-diagonal terms; in par- 
ticular, the off-diagonal term 



Cov 



G-GWT^-G 



2GY 



: (C G-G + Gf) 



(21 + l)/ sky 



(7) 



is non-zero. Inclusion of this term is necessary to avoid 
"double-counting" information in the final analysis (a re- 
lated subtlety is that the filter defined in Eq. [S] should 
not be varied from the fiducial value when computing 
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parameter 


fiducial value 


A error 


B error 




0.1400 


0.0016 


0.0016 


fl b h 2 


0.02400 


0.00019 


0.00019 


^LSS 


1.390 Gpc 


0.029 Gpc 


0.023 Gpc 




0.84 


0.10 


0.054 


w 


-1.000 


0.099 


0.081 


dw/da 


0.00 


0.18 


0.18 


"c" 


0.00 


1.1 


0.48 




0.350 Mpc" 1 


0.071 Mpc" 1 


0.062 Mpc' 1 




17.00 


0.18 


0.17 


b 


1.00 


0.18 


0.10 



Fig. 3. — Computed errors from a Fisher matrix analysis combin- 
ing WMAP, ACT and SALT datasets as described. Set "A" con- 
tains the temperature anisotropy and galaxy count signals, set "B" 
includes the cross-correlation between the temperature squared and 
galaxy count signals. 



the derivatives dCi/dX.) An additional term in Eq. [7| 
proportional to the three-point correlation of the galaxy 
signal multiplied by the total temperature RMS, we have 
neglected as being strongly subdominant. 

Other off-diagonal terms are zero, at least in the (good) 
approximation that the (unsquared) temperature-galaxy 
correlation is zero; we expect the lensing signal to be 
much smaller than the kSZ. Since the SALT galaxy sam- 
ple only covers 100 deg 2 , we can ignore the (very) large- 
scale correlations induced by the integral Sachs- Wolfe ef- 
fect which were detected in a cross-correlation of WMAP 
and SDSS data. Since the covariance matrix is block di- 
agonal, its inversion to produce the Fisher matrix is then 
trivial. 

7. ALL-SCALE LIMITS 

To demonstrate for which parameters a kSZ-galaxy 
cross-correlation provides the greatest information, we 
give the errors for two separate analyses, one that uses 
WMAP and ACT temperature anisotropy, and SALT 
galaxy count, data (set "A"), and one that uses this data 
plus the cross-correlation signal between the ACT and 
SALT data (set "B"). 

We first consider constraints that can be made by in- 
cluding information on all scales - i.e., when we use all 
available I measurements limited only by Poisson and 
instrument noise. Our results are listed in the table in 
Fig. 3, which gives the fiducial cosmological model and 
estimated errors on our suite of parameters. 

For simplicity, we have assumed a flat universe prior. 
Because of recent interest in dark energy models that 
display "phantom" behaviour, we do not require w > — 1; 
however, we assume that the dark energy fluid does not 
lead to appreciable clustering, as is the case, e.g., when 
the true sound speed c s is unity. 

As can be seen, the addition of the cross-correlation 
information reduces the error on as, cZlss and w, and 
gives significant constraints on the parameters associated 
with late-time gas physics, c and cr sm . 

The sources of these results are not hard to understand. 
The strong erg dependence i n the cross-correla tion signal 
has already been noted bv iDore et alJ l)2004|) . and the 
reduction in error on c?lss and w are presumably due to 
the dependence of the velocity field on the instantaneous 
value of the Hubble parameter, Eq. ^ 

That errors on dw/da are not significantly improved 



is most likely due to our choice of parametrization; a 
dark energy model that restricted variations in w to very 
late times (z < 1) would presumably receive much bet- 
ter constraints from the addition of the cross-correlation 
information. The redshift of reionization, z r i, is strongly 
constrained by the overa ll amplitude of the kSZ signal, 
as noted bv lZhang et al.1 l)2004|) ; a slight improvement in 
limits on including the cross-correlation comes from re- 
moving degeneracies with other parameters, in particular 
its- Meanwhile, the significant improvement in parame- 
ters c and <j sm is unsurprising; cross-correlation is an 
excellent method to study the late-time time evolution 
and spatial dependence of the ionized baryons. 

8. LINEAR VS. NON-LINEAR SCALES 

While we have attempted to address a collection of 
different effects that complicate the analysis of galaxy 
and kSZ information, there remain a number of diffi- 
culties that will complicate any future analysis. A sig- 
nificant and troublesome systematic effect, that comes 
in at small scales, is the emergence of non-linear, and 
poss ibly non-deterministic, bias in the galaxy count sig- 
nal ijScherrer fc Weinberdll9 98). Our analysis above ne- 
glected these effects. 

As can be seen in Eq. ^ the velocity signal is pro- 
portional to 5/k~ 2 , and thus is stronger at larger scales 
where bias models are better constrained. It can thus 
be said that the kSZ signal, in addition to having differ- 
ent systematics from a galaxy-count signal, also provides 
"higher-quality" information that may be less susceptible 
to the difficulties inherent in understanding structures on 
very small scales. 

We can quantify this suggestion by redoing the anal- 
ysis in the previous section and discarding small-scale 
information. Since the bulk of our galaxies are found at 
z « 0.7, this can be done simply by including only low-£ 
information. We thus consider a cut at £ ss 1750, corre- 
sponding to a comoving scale of lh Mpc _1 at z of 0.7. In 
general, effects such as a scale dependence in the bias are 
expected to become strong at this scale; it is important 
to note that such effects do appear to a lesser extent at 
even larger scales ■' ISlanton et al.lll99^ . 

The table in Fig. 4 shows the results of this cut. As can 
be seen, when the signal is limited to larger scales, the 
kSZ signal indeed becomes more important in constrain- 
ing cosmological parameters. A "conservative" analy- 
sis of the WMAP, ACT and SALT data that discards 
smaller scales will improve its errors by a factor of two or 
more when it includes information derived from a cross- 
correlation of the kSZ and galaxy count signals. 

As mentioned above, the bias is not expected to be 
completely scale-free linear at any of the scales probed 
by an ACT/SALT experiment. Furthermore, even when 
we restrict our study to k < lh Mpc -1 , we are still in the 
"non-linear" regime in one important sense; below scales 
of O.lh Mpc -1 (corresponding to I of 175), the dark mat- 
ter power spectrum can no longer be well-approximated 
by linear theory alone. Any analysis of either the kSZ or 
galaxy signal must always take this kind of non-linearity 
into account. 

9. DISCUSSION 

One of the lessons of our paper is that the kSZ is of 
interest not only for studies of dark energy, but also 
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parameter 


Ih Mpc 1 A error 


Ih Mpc 1 B error 


n m ti 2 


0.0021 


0.0017 


n b h 2 


0.00043 


0.00020 




0.032 Gpc 


0.029 Gpc 


°"8 


0.55 


0.22 


W 


0.11 


0.10 


dw/da 


0.19 


0.19 


"c" 


n.c. 


1.4 




0.81 Mpc" 1 


0.14 Mpc" 1 




0.39 


0.19 


b 


0.72 


0.31 



Fig. 4. — The effect of restricting the signal to large scales 
where the bias is more easily modeled. We truncate the signal 
at lh Mpc -1 and consider both the case in which the CMB and 
galaxy-count signals are used only separately ("set A"), and the 
"set B" analysis, where the cross-correlation information is in- 
cluded. The term "n.c." is used when constraints are so loose as to 
signal a significantly non-Gaussian distribution where parameters 
refer to physical quantities that may not go negative. 



In particular, iHernandez-Monteagudo et alJ ([2005) 
uses the kSZ to provide an unbiased measurement of in- 
dividual cluster motions, whereas we consider the pos- 
sibility of a statistical extraction of velocity fields on a 
wide range of scales. The physics behind the generation 
of peculiar velocities is rich, and a number of different 
kinds of analysis made be done to produce results with 
different systematics. 

Both a galaxy-count and a galaxy-kSZ cross- 
correlation have excellent prospects for the measurement 
of cosmological parameters, but have different systemat- 
ics and are sensitive to different physics. As we have 
discussed in Sec. |H1 the kSZ can be said to provide "high 
quality" information, by relying for the bulk of its signal 
on velocity flows on large scales where a cross-correlation 
does not have to worry as much about the effects of non- 
linear and non-deterministic bias. 



for investigations of the distribution and evolution of 
baryons in the universe. Our work is comp lementary 
to that of lHernandez-Monteagudo et al.1 l|2005j) : both pa- 
pers consider the possibilities of using the ACT measure- 
ment of the kSZ to study cosmological parameters, but 
do so in different ways. 
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